A blood test to monitor bee health across a European network of agricultural sites of different land-use by MALDI BeeTyping mass spectrometry.

There are substantial concerns about impaired honey bee health and colony losses due to several poorly understood factors. We used MALDI profiling (MALDI BeeTyping®) analysis to investigate how some environmental and management factors under field conditions across Europe affected the honey bee haemolymph peptidome (all peptides in the circulatory fluid), as a profile of molecular markers representing the immune status of Apis mellifera. Honey bees were exposed to a range of environmental stressors in 128 agricultural sites across eight European countries in four biogeographic zones, with each country contributing eight sites each for two different cropping systems: oilseed rape (OSR) and apple (APP). The full haemolymph peptide profiles, including the presence and levels of three key immunity markers, namely the antimicrobial peptides (AMPs) Apidaecin, Abaecin and Defensin-1, allowed the honey bee responses to environmental variables to be discriminated by country, crop type and site. When considering just the AMPs, it was not possible to distinguish between countries by the prevalence of each AMP in the samples. However, it was possible to discriminate between countries on the amounts of the AMPs, with the Swedish samples in particular expressing high amounts of all AMPs. A machine learning model was developed to discriminate the haemolymphs of bees from APP and OSR sites. The model was 90.6 % accurate in identifying the crop type from the samples used to build the model. Overall, MALDI BeeTyping® of bee haemolymph represents a promising and cost-effective "blood test" for simultaneously monitoring dozens of peptide markers affected by environmental stressors at the landscape scale, thus providing policymakers with new diagnostic and regulatory tools for monitoring bee health.

Novel indices reveal that pollinator exposure to pesticides varies across biological compartments and crop surroundings.

Declines in insect pollinators have been linked to a range of causative factors such as disease, loss of habitats, the quality and availability of food, and exposure to pesticides. Here, we analysed an extensive dataset generated from pesticide screening of foraging insects, pollen-nectar stores/beebread, pollen and ingested nectar across three species of bees collected at 128 European sites set in two types of crop. In this paper, we aimed to (i) derive a new index to summarise key aspects of complex pesticide exposure data and (ii) understand the links between pesticide exposures depicted by the different matrices, bee species and apple orchards versus oilseed rape crops. We found that summary indices were highly correlated with the number of pesticides detected in the related matrix but not with which pesticides were present. Matrices collected from apple orchards generally contained a higher number of pesticides (7.6 pesticides per site) than matrices from sites collected from oilseed rape crops (3.5 pesticides), with fungicides being highly represented in apple crops. A greater number of pesticides were found in pollen-nectar stores/beebread and pollen matrices compared with nectar and bee body matrices. Our results show that for a complete assessment of pollinator pesticide exposure, it is necessary to consider several different exposure routes and multiple species of bees across different agricultural systems.

Distribution of infectious and parasitic agents among three sentinel bee species across European agricultural landscapes.

Infectious and parasitic agents (IPAs) and their associated diseases are major environmental stressors that jeopardize bee health, both alone and in interaction with other stressors. Their impact on pollinator communities can be assessed by studying multiple sentinel bee species. Here, we analysed the field exposure of three sentinel managed bee species (Apis mellifera, Bombus terrestris and Osmia bicornis) to 11 IPAs (six RNA viruses, two bacteria, three microsporidia). The sentinel bees were deployed at 128 sites in eight European countries adjacent to either oilseed rape fields or apple orchards during crop bloom. Adult bees of each species were sampled before their placement and after crop bloom. The IPAs were detected and quantified using a harmonised, high-throughput and semi-automatized qPCR workflow. We describe differences among bee species in IPA profiles (richness, diversity, detection frequencies, loads and their change upon field exposure, and exposure risk), with no clear patterns related to the country or focal crop. Our results suggest that the most frequent IPAs in adult bees are more appropriate for assessing the bees' IPA exposure risk. We also report positive correlations of IPA loads supporting the potential IPA transmission among sentinels, suggesting careful consideration should be taken when introducing managed pollinators in ecologically sensitive environments.

Pesticide use negatively affects bumble bees across European landscapes.

Sustainable agriculture requires balancing crop yields with the effects of pesticides on non-target organisms, such as bees and other crop pollinators. Field studies demonstrated that agricultural use of neonicotinoid insecticides can negatively affect wild bee species1,2, leading to restrictions on these compounds3. However, besides neonicotinoids, field-based evidence of the effects of landscape pesticide exposure on wild bees is lacking. Bees encounter many pesticides in agricultural landscapes4-9 and the effects of this landscape exposure on colony growth and development of any bee species remains unknown. Here we show that the many pesticides found in bumble bee-collected pollen are associated with reduced colony performance during crop bloom, especially in simplified landscapes with intensive agricultural practices. Our results from 316 Bombus terrestris colonies at 106 agricultural sites across eight European countries confirm that the regulatory system fails to sufficiently prevent pesticide-related impacts on non-target organisms, even for a eusocial pollinator species in which colony size may buffer against such impacts10,11. These findings support the need for postapproval monitoring of both pesticide exposure and effects to confirm that the regulatory process is sufficiently protective in limiting the collateral environmental damage of agricultural pesticide use.

Sulfoxaflor effects depend on the interaction with other pesticides and Nosema ceranae infection in the honey bee (Apis mellifera).

Honey bees health is compromised by many factors such as the use of agrochemicals in agriculture and the various diseases that can affect them. Multiple studies have shown that these factors can interact, producing a synergistic effect that can compromise the viability of honey bees. This study analyses the interactions between different pesticides and the microsporidium Nosema ceranae and their effect on immune and detoxification gene expression, sugar consumption and mortality in the Iberian western honey bee (Apis mellifera iberiensis). For this purpose, workers were infected with N. ceranae and subjected to a sugar-water diet with field concentrations of the pesticides sulfoxaflor, azoxystrobin and glyphosate. Increased sugar intake and altered immune and cytochrome P450 gene expression were observed in workers exposed to sulfoxaflor and infected with N. ceranae. None of the pesticides affected Nosema spore production in honey bee gut. Of the three pesticides tested (alone or in combination) only sulfoxaflor increased mortality in honey bees. Taken together, our results suggest that the effects of sulfoxaflor were attenuated in contact with other pesticides, and that Nosema infection leads to increase sugar intake in sulfoxaflor-exposed bees. Overall, this underlines the importance of studying the interaction between different stressors to understand their overall impact not only on honey bee but also on wild bees health.

Evolutionarily diverse origins of deformed wing viruses in western honey bees.

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

"Migratory beekeeping and its influence on the prevalence and dispersal of pathogens to managed and wild bees".

Demand for food is growing along with the human population, leading to an increase in plant production. Many crops are pollinated by insects, so the global demand for managed pollinators is also increasing. The honey bee has traditionally been considered the main provider of crop pollination services. For providing it beekeepers seasonally transport hives to different locations after the flowering of different crops. These movements could be detrimental to pollinators by: i) stressing honey bees, making them more susceptible to pathogens and parasites; ii) spreading bee parasites and pathogens across locations; iii) increasing the transmission of parasites and pathogens between managed and wild pollinators and vice versa (spillover and spillback, respectively). To understand the impact of migratory beekeeping on bee health, we conducted a systematic review to identify the main trends and provide a complete picture of existing knowledge on the subject. We found 52 studies analysing pathogen-related impacts of migratory beekeeping on honey bees. However, only 16 investigations tested the effect of migratory practices on the prevalence and spread of pathogens and parasites. We found no studies that assessed the impact of migratory beekeeping on the occurrence and spread of pests and diseases in wild bees. In general, migratory beekeeping tends to increase the prevalence of pathogens and parasites in honey bee colonies. However, the results were very heterogeneous, probably due to several uncontrolled underlying factors such as management, biological and geographical factors, and the interactions between them. In conclusion, there is an urgent need for studies to assess the impact of migratory beekeeping on bee health, given the current global bee decline and the expected increase in migratory beekeeping due to climate change and crop pollination demand.

Honey bees and climate explain viral prevalence in wild bee communities on a continental scale.

Viruses are omnipresent, yet the knowledge on drivers of viral prevalence in wild host populations is often limited. Biotic factors, such as sympatric managed host species, as well as abiotic factors, such as climatic variables, are likely to impact viral prevalence. Managed and wild bees, which harbor several multi-host viruses with a mostly fecal-oral between-species transmission route, provide an excellent system with which to test for the impact of biotic and abiotic factors on viral prevalence in wild host populations. Here we show on a continental scale that the prevalence of three broad host viruses: the AKI-complex (Acute bee paralysis virus, Kashmir bee virus and Israeli acute paralysis virus), Deformed wing virus, and Slow bee paralysis virus in wild bee populations (bumble bees and solitary bees) is positively related to viral prevalence of sympatric honey bees as well as being impacted by climatic variables. The former highlights the need for good beekeeping practices, including Varroa destructor management to reduce honey bee viral infection and hive placement. Furthermore, we found that viral prevalence in wild bees is at its lowest at the extreme ends of both temperature and precipitation ranges. Under predicted climate change, the frequency of extremes in precipitation and temperature will continue to increase and may hence impact viral prevalence in wild bee communities.

Detection of Microsporidia in Pollinator Communities of a Mediterranean Biodiversity Hotspot for Wild Bees.

Insect pollination is crucial for the maintenance of natural and managed ecosystems but the functioning of this ecosystem service is threatened by a worldwide decline of pollinators. Key factors in this situation include the spread and interspecific transmission of pathogens worldwide through the movement of managed pollinators. Research on this field has been mainly conducted in some particular species, while studies assessing the interspecific transmission of pathogens at a community level are scarce. However, this information is pivotal to design strategies to protect pollinators. Herein, we analysed the prevalence of two common microsporidia pathogens of managed honey bees (Nosema ceranae and N. apis) in bee communities of semiarid Mediterranean areas from the Southeast of the Iberian Peninsula. Our results confirm the ability of N. ceranae to disperse across wild bee communities in semiarid Mediterranean ecosystems since it was detected in 36 Apoidea species (39% of the sampling; for the first time in nine genera). The prevalence of the pathogen did not show any phylogenetic signal which suggests a superfamily host range of the pathogen or that wild bees may be acting only as vectors of N. ceranae. In addition, N. apis was detected in an Eucera species, which is the second time it has been detected by molecular techniques in a host other than the honey bee. Our study represents the primary assessment of the prevalence of microsporidia at community level in Mediterranean areas and provides outstanding results on the ability of Nosema pathogens to spread across the landscape.

Bee Trypanosomatids: First Steps in the Analysis of the Genetic Variation and Population Structure of Lotmaria passim, Crithidia bombi and Crithidia mellificae.

Trypanosomatids are among the most prevalent parasites in bees but, despite the fact that their impact on the colonies can be quite important and that their infectivity may potentially depend on their genotypes, little is known about the population diversity of these pathogens. Here we cloned and sequenced three non-repetitive single copy loci (DNA topoisomerase II, glyceraldehyde-3-phosphate dehydrogenase and RNA polymerase II large subunit, RPB1) to produce new genetic data from Crithidia bombi, C. mellificae and Lotmaria passim isolated from honeybees and bumblebees. These were analysed by applying population genetic tools in order to quantify and compare their variability within and between species, and to obtain information on their demography and population structure. The general pattern for the three species was that (1) they were subject to the action of purifying selection on nonsynonymous variants, (2) the levels of within species diversity were similar irrespective of the host, (3) there was evidence of recombination among haplotypes and (4) they showed no haplotype structuring according to the host. C. bombi exhibited the lowest levels of synonymous variation (πS= 0.06 ± 0.04 %) - and a mutation frequency distribution compatible with a population expansion after a bottleneck - that contrasted with the extensive polymorphism displayed by C. mellificae (πS= 2.24 ± 1.00 %), which likely has a more ancient origin. L. passim showed intermediate values (πS= 0.40 ± 0.28 %) and an excess of variants a low frequencies probably linked to the spread of this species to new geographical areas.

Spatial and temporal patterns of genetic diversity in Bombus terrestris populations of the Iberian Peninsula and their conservation implications.

The bumblebee Bombus terrestris is used worldwide for crop pollination. Despite its positive impact on crop yield, it has become a widespread threat to biodiversity due to its interactions with local bumblebee populations. Commercial subspecies introduced to the Iberian Peninsula since the 1990s without any regulation have colonized the environment, with evidence of naturalization and introgression with the endemic subspecies Bombus terrestris lusitanicus. We have used mitochondrial and nuclear genetic data to describe the current genetic diversity of the Iberian population and to estimate the expansion of commercial bumblebees. Samples from the natural distribution range of the commercial subspecies, the natural intergradation area between the two subspecies and from a period prior to the use of commercial colonies (i.e., before the 1990s) have been used for comparison. Our results show that the mitochondrial haplotype of the commercial breeds has spread throughout the territory, which, together with subtle changes observed in the nuclear genetic diversity of the populations, indicates that hybridization and consequent introgression are occurring in most of the peninsula. It is, therefore, necessary to improve the existing legislation concerning the management and exportation of commercial bumblebees to conserve locally adapted populations.

Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs.

BACKGROUND: With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and FST) to select the most informative SNPs for ancestry inference. RESULTS: Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof. CONCLUSIONS: The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

Wide diversity of parasites in Bombus terrestris (Linnaeus, 1758) revealed by a high-throughput sequencing approach.

Assessing the extent of parasite diversity requires the application of appropriate molecular tools, especially given the growing evidence of multiple parasite co-occurrence. Here, we compared the performance of a next-generation sequencing technology (Ion PGM ™ System) in 12 Bombus terrestris specimens that were PCR-identified as positive for trypanosomatids (Leishmaniinae) in a previous study. These bumblebees were also screened for the occurrence of Nosematidae and Neogregarinorida parasites using both classical protocols (either specific PCR amplification or amplification with broad-range primers plus Sanger sequencing) and Ion PGM sequencing. The latter revealed higher parasite diversity within individuals, especially among Leishmaniinae (which were present as a combination of Lotmaria passim, Crithidia mellificae and Crithidia bombi), and the occurrence of taxa never reported in these hosts: Crithidia acanthocephali and a novel neogregarinorida species. Furthermore, the complementary results produced by the different sets of primers highlighted the convenience of using multiple markers to minimize the chance of some target organisms going unnoticed. Altogether, the deep sequencing methodology offered a more comprehensive way to investigate parasite diversity than the usual identification methods and provided new insights whose importance for bumblebee health should be further analysed.

The Effect of Migratory Beekeeping on the Infestation Rate of Parasites in Honey Bee (Apis mellifera) Colonies and on Their Genetic Variability.

Migratory beekeeping is a widely extended practice aimed at increasing the yield of products and pollination services of honey bee colonies. However, it represents a stress factor, as it facilitates the dissemination of diseases and may compromise the genetic identity of the colonies involved. To analyze the extent of these effects, pathogens infestation rate and genetic composition were monitored in a field experiment comparing stationary and migratory colonies sharing the same environmental conditions but differing in management (stationary vs. migratory) and genetic background. We studied the pathogens infestation rate (Varroa destructor, Nosema spp., and Deformed Wing Virus (DWV)) at four different times: before migratory operation, two weeks later, at the end of the migratory period, and two weeks after the return of the migratory hives. An increased incidence of V. destructor and Nosema ceranae and a lower DWV viral load were found in migratory colonies. Temporary changes in genetic diversity were detected regardless of colony type, suggesting that stressors other than management affect the genetic diversity of the colonies. Our study demonstrates that migratory practices have variable effects on the health and genetic diversity of honey bee colonies, which should be taken into account for the development of sustainable beekeeping.

Intercontinental long-distance seed dispersal across the Mediterranean Basin explains population genetic structure of a bird-dispersed shrub.

Long-distance dispersal (LDD) is a pivotal process for plants determining their range of distribution and promoting gene flow among distant populations. Most fleshy-fruited species rely on frugivorous vertebrates to disperse their seeds across the landscape. While LDD events are difficult to record, a few ecological studies have shown that birds move a sizeable number of ingested seeds across geographic barriers, such as sea straits. The foraging movements of migrant frugivores across distant populations, including those separated by geographic barriers, creates a constant flow of propagules that in turn shapes the spatial distributions of the genetic variation in populations. Here, we have analysed the genetic diversity and structure of 74 populations of Pistacia lentiscus, a fleshy-fruited shrub widely distributed in the Mediterranean Basin, to elucidate whether the Mediterranean Sea acts as a geographic barrier or alternatively whether migratory frugivorous birds promote gene flow among populations located on both sides of the sea. Our results show reduced genetic distances among populations, including intercontinental populations, and they show a significant genetic structure across an eastern-western axis. These findings are consistent with known bird migratory routes that connect the European and African continents following a north-southwards direction during the fruiting season of many fleshy-fruited plants. Further, approximate Bayesian analysis failed to explain the observed patterns as a result of historical population migrations at the end of Last Glacial Maximum. Therefore, anthropic and/or climatic changes that would disrupt the migratory routes of frugivorous birds might have genetic consequences for the plant species they feed upon.

The toxic unit approach as a risk indicator in honey bees surveillance programmes: A case of study in Apis mellifera iberiensis.

The influence of genetic diversity and exposure to xenobiotics on the prevalence of pathogens was studied within the context of a voluntary epidemiological study in Spanish apiaries of Apis mellifera iberiensis, carried out during the spring season of years 2014 and 2015. As such, the evolutionary lineages of the honey bee colonies were identified, a multiresidue analysis of xenobiotics was carried out in beebread and worker bee samples, and the Toxic Unit (TUm) was estimated for each sampled apiary. The relationship between lineages and the most prevalent pathogens (Nosema ceranae, Varroa destructor, trypanosomatids, Black Queen Cell Virus; and Deformed Wing Virus) was analysed with contingency tables, and the possible relationships between TUm and the prevalence of these pathogens were studied by using a factor analysis. The statistical analysis supported the associations between V. destructor and Deformed Wing Virus (DWV), and between N. ceranae and Black Queen Cell Virus (BQCV), but the association between these pathogens and trypanosomatids was not observed. TUm values varied between 5.5 × 10-6 and 3.65 × 10-1. When TUm < 3.35 × 10-4, it was mainly determined by coumaphos, tau-fluvalinate and/or chlorfenvinphos. At higher values, other insecticides also contributed to TUm, although a clear predominance was not seen up to TUm ≥ 1.83 × 10-2, when it was mainly defined by acrinathrin, spinosad and/or imidacloprid. The possible cumulative effect from the joint action of xenobiotics was >10% in the 63% of the cases. The prevalence of pathogens did not appear to be influenced by the distribution of evolutionary lineages and, while the prevalence of V. destructor was not found to be determined by TUm, there was a trend towards an increasing prevalence of N. ceranae when TUm ≥ 23 10-4. This study is an example of using TUm approach beyond the field of the ecotoxicology.

An integrative approach to discern the seed dispersal role of frugivorous guilds in a Mediterranean semiarid priority habitat.

Seed dispersal is an essential process to maintain the viability of plant populations, and understanding this ecological process allows management strategies to be developed to conserve ecosystems. European Union priority habitat 5220* is defined as "Mediterranean arborescent shrubland with Ziziphus lotus" and it represents a favorable microclimate within the severe climatic conditions typical of the semiarid south-eastern region of the Iberian Peninsula. Therefore, the study of seed dispersal in this priority habitat by different frugivorous guilds, is a challenge for its conservation. In this study, we have characterized a mutualistic network of seed dispersal that is mediated by vertebrates (mammals and birds) in the protected habitat 5220*. The aims of this study were to: (i) identify the seed disperser community; (ii) analyze the relative role of key species in the dispersal process; and (iii) compare the functional ecology of the seed dispersal process between mammals and birds. As such, we collected animal faeces to determine seed dispersers taxonomy, identifying the mammals through the visual aspect of the faeces and the birds by DNA barcoding. In the case of birds, we also collected regurgitated seeds in which the disperser species was also identified through molecular techniques. This allowed us to build-up a mutualistic network and to identify the relative role of these animals in seed dispersal. Our results showed that mammals and birds fulfilled complementary roles in seed dispersal, with birds representing the main dispersers of key plants within the 5220* habitat, and mammals the main dispersers of human-cultivated plants. Herein, we provide a useful approach with relevant information that can be used to propose management policies that focus on restoring the threatened 5220* habitat, promoting the role of birds to disperse key species that structure plant communities of this priority habitat.

Coexistence of genetically different Varroa destructor in Apis mellifera colonies.

The aim of this study was to investigate the genetic diversity of Varroa destructor parasitizing Apis mellifera colonies and to test for possible host-parasite association at the mitochondrial DNA (mtDNA) level. Six A. mellifera haplotypes (including a novel C2aa) and five haplotypes of V. destructor were detected in 29 analyzed colonies from eight sampling sites in Serbia. We revealed the presence of the K and S1 haplotypes as well as KS1 and KP1 heteroplasmic mite individuals in all localities, while the P1 haplotype was only found in four sampling sites. Significant differences in V. destructor genetic diversity were found at both apiary and colony levels, with mite haplotypes coexisting in almost all tested colonies. In addition, a significant correlation between the number of analyzed mites per colony and the number of identified V. destructor haplotypes was observed. However, no significant host-parasite relationship was found, suggesting that mites bearing different haplotypes as well as those heteroplasmic individuals are well adapted to the host, A. mellifera, independently of the identified haplotype present in each colony. Our results will contribute to future population and biogeographic studies concerning V. destructor infesting A. mellifera, as well as to better understanding their host-parasite relationship.

A new multiplex PCR protocol to detect mixed trypanosomatid infections in species of Apis and Bombus.

Trypanosomatids are highly prevalent pathogens of Hymenoptera; however, most molecular methods used to detect them in Apis and Bombus spp. do not allow the identification of the infecting species, which then becomes expensive and time consuming. To overcome this drawback, we developed a multiplex PCR protocol to readily identify in a single reaction the main trypanosomatids present in these hymenopterans (Lotmaria passim, Crithidia mellificae and Crithidia bombi), which will facilitate the study of their epidemiology and transmission dynamics. A battery of primers, designed to simultaneously amplify fragments of the RNA polymerase II large subunit (RPB1) of L. passim, the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of C. mellificae and the DNA topoisomerase II (TOPII) of C. bombi, was tested for target specificity under single and mixed template conditions using DNA extracted from cell cultures (L. passim ATCC PRA403; C. mellificae ATCC 30254) and from a bumblebee specimen infected with C. bombi only (14_349). Once validated, the performance of the method was assessed using DNA extractions from seven Apis mellifera (Linnaeus, 1758) and five Bombus terrestris (Linnaeus, 1758) field samples infected with trypanosomatids whose identity had been previously determined by PCR-cloning and sequencing (P-C-S). The new method confirmed the results obtained by P-C-S: two of the honeybee samples were parasitized by L. passim, C. mellificae and C. bombi at the same time, whereas the other five were infected with L. passim only. The method confirmed the simultaneous presence of L. passim and C. mellificae in two B. terrestris, where these parasites had not previously been reported.

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera).

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M-lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C-lineage). Effective conservation actions require reliable molecular tools to identify pure-bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q-values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.